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ABSTRACT 

The 6 X 10® Mq supermassive black hole at the center of the giant elliptical galaxy M87 powers a relativistic 
jet. Observations at millimeter wavelengths with the Event Horizon Telescope have localized the emission from 
the base of this jet to angular scales comparable to the putative black hole horizon. The jet might be powered 
directly by an accretion disk or by electromagnetic extraction of the rotational energy of the black hole. However, 
even the latter mechanism requires a confining thick accretion disk to maintain the required magnetic flux near the 
black hole. Therefore, regardless of the jet mechanism, the observed jet power in M87 implies a certain minimum 
mass accretion rate. If the central compact object in M87 were not a black hole but had a surface, this accretion 
would result in considerable thermal near-infrared and optical emission from the surface. Current flux limits on 
the nucleus of M87 strongly constrain any such surface emission. This rules out the presence of a surface and 
thereby provides indirect evidence for an event horizon. 

Subject headings: black hole physics - galaxies: individual (M87) - gravitation - radio continuum: galaxies - 
infrared: galaxies - ultraviolet: galaxies 


1. INTRODUCTION 

It is now widely accepted that active galactic nuclei (AGN) 
are powered by supermassive objects (reaching 10 ^®Mq) that 
are sufficiently compact to exclud e all other astrophysically 
credible alternatives to black holes (lReeslll9^ . However, it is 
less clear that these objects possess the defining characteristic 
of a black hole: an event horizorQ. The existence of black 
hole event horizons plays a central role in a number of 
puzzles associated with black holes, e.g., the information 
paradox. A number of recent results suggest that a resolution 
of th ese puzzles may result in modifications on horizo n scales 
(e.g.. lMathudl201 lllAlmheiri et al.llT013HMathu^l2014^ ■ which 
provides strong motivation for seeking astronomical evidence 
for the reality of event horizons. 

Accretion onto compact objects with a surface, e.g., white 
dwarfs, neutrons stars, results in the formation of a boundary 
layer in which any remaining kinetic energy contained within 
the accretion flow is thermalized and radiated. In contrast, gas 
accreting onto a black hole is free to advect any excess energy 
across the horizon without further observational consequence. 
If the mass accretion rate, M, can be independently estimated, 
this difference provides an observational means to distinguish 
between the presence of a surface, or more accurately a 
“photosphere,” and a horizon. 

The above argument has already been used to argue 

* Here we will employ an astrophysically motivated definition of the 
horizon: a surface from inside which astronomical signals cannot propagate 
to large distances in astronomically relevant timescales. Formally, for a 
dynamical system, such a surface is identified with the apparent horizon. 
However, in the context of astrophysical black holes described by general 
relativity, it corresponds to the event horizon as well. 


for the existence of event horizons in X-ray binaries by 
com paring neutron star and black hole systems in aggre- 


gate (iNaravan et al.lll997 

;lGarcia et al.ll200lHNaravan & Hey 11 

120021 iDone & Gierlihskil 

2003HNaravan & Mcnintockll2008li. 


However, the advent of horizon-resolving observations, en¬ 
abled by millimeter-wavelength very long baseline interero- 
metric observations ( mm-VLBI) carried out by the Event Hori- 


zon Telescope (EHT. iDoeleman et al. 

120091 lDoe1emanll2010l 

IDoeleman et al.ll2008t lEish et alJI201 1 

IDoeleman et alJI2012l). 


has made it possible to extended the argument to individual 
systems. This is primarily because restricting the size of any 
photospheric emission to horizon scales enables two important 
simplifications: 

1. Any putative photosphere that lies within the photon 
orbit is expected to radiate to a good approximation 
like a black body, independently of the details of its 
composition (iBroderick & Naravanll200^ 1200711 . This 
is because a majority of the photons emitted from 
the photosphere will be strongly lensed back onto the 
photosphere, thermally coupling the photosphere to 
itself and to the emitted photon held. As the redshift 
of the surface increases, the blackbody approximation 
becomes increasingly accurate. 

2. The expected temperature of the photosphere emission, 
as seen by distant observers, is dependent upon the 
mass accretion rate M and the apparent photosphere size, 
the latter of which is fixed when the photosphere lies 
within the photon orbit. Thus, assuming that the system 
has reached steady stat^E any independent estimate of 

^ The additional gravitational time delay for radiation to escape from 





















2 


BRODERICK ET AL. 


M immediately determines both the luminosity and the 
radiation spectrum as seen by a distant observer. 

Essentially, by restricting the surface to be sufficiently 
compact, it is possible to robustly predict the properties of 
any putative photosphere emission, independent of the specific 
properties of the radiating surface. Direct flux limits can 
then be used to constrain and/or exclude the presence of a 
photosphere. 

The above argument has already been successfully employed 
in the case of Sagittarius A* (Sgr A*), th e supermassive black 
hole at the center of the M i lky W a y dNaravan et al.l 119981: 
iBroderick & NarayMil 120061 120071 : iNaravan & McClintocij 
120081 : iBroderick et'^l2009l) . where limits on M were obtained 
from the observed bolometric luminosity, assumed to arise 
from the accretion flow during infall. It was shown that, for 
physically reasonable radiative efficiencies, it is impossible to 
accommodate a photosphere. Therefore, Sgr A* must have an 
event horizon behind which the kinetic energy of the infalling 
accreting gas is hidden. 

The supermassive black hole at the center of the nearby 
giant elliptical galaxy M87 (which we will refer to as M87* 
hereafter) is a second target for which the EHT has provided 
horizon-scale limits upon the extent of its mm-wavelength 
emission. Here we explore the implications of these and 
related observations for the existence of an event horizon in 
M87*. In what follows we assume a mass of 6.16 x 10 ®Mq 
and dista nce of 16.5 Mpc for M 87*, reported in the recent 
review bv iKormendv & (120131) . The mass e stimate is based 

on the stellar dynamical modeling described in iGebhardt et alJ 
dMl, and is roughly a factor of two larg er than the 
value obtained by gas dynamical measurements (iWalsh et al.l 
120131) : both methods are potentially complicated by the fact 
that M87* is offset from the center-of-light by as much as 
10 pc, im pacting the underlying a ssumptions regarding orbital 
isotropy (IBatcheldor et al.l 120101) . However, our qualitative 
conclusions are insensitive to this difference, being marginally 
stronger with the smaller black hole mass0 

M87* is nearly three orders of magnitude more massive than 
Sgr A*, probing a mass regime more relevant for the bright 
AGN observed at high redshift. Unlike Sgr A*, whose mm and 
radio emission seem to be primarily from the accretion flow 
dYuan & Naravanll2014ft . the mm/radio emission of M87* is 
dominated by a powerful relativistic jet. 

Two lines of argument strongly imply that the relativistic 
jet in M87 originates near the black hole. Eirst, astrometric 
measurements of the radio core position reveal a wavelength- 
dependent shift, asymptoting to a fixed position at short wave¬ 
length s, consistent with a black-hole-launched jet (iHada et al.l 
1201 ll) . Second, the small scales implied by EHT observations 
at 1.3mm are commensurate with the scales expected near the 
black hole (iDoeleman et al.l2012L and references therein). This 
is additionally supported by the success of semi-analytical jet 
models in simultaneously reproducing the EHT observations, 
large-scale jet properties, and spectral energy density (SED) of 
M87, that will be reported elsewhere. Thus, the sub-mm jet is 

a compact surface is insufficient to prevent the system from reaching 
steady state, since the timescale diverges only logarithmically as the radius 
of the surface approaches the horizon radius IBroderick & NaravanI I20Q6I: 
INaravan & McClintockll20Q^ : IBroderick et alJ2009n . 

^ The temperature of a putative surface observed at infinity, defined in 
Equation ED depends on tlie mass as M ^^d tFius is Fiigher for smaller 
masses. This shifts the resulting emission to higher frequencies where better 
limits on the observed nuclear emission exist. 


almost certainly launched in the immediate environment of the 
black hole. 

All current potential mechanisms for launching relativistic 
jets invoke an accretion flo w. This is trivially tru e for 
disk-launched outflows (e.g.. iBlandford & Pavnelfl982h . but 
it is also the case for b lack hole spin-powered jets (e.g., 
IBlandford & Znai^ll977h . In the latter case, the accretion 
disk is needed to confine the horizon-penetrating magnetic 
flux which enables the black hole’s rotational energy to be 
tapped. While jet efficiencies, defined as ryjet = L^ei/Mc^ 
where Ljet is the total jet luminosity (radiative, magnetic 
and mechanical), can instantaneously exceed unity as a 
result of the electromagnetic tapping of the black hole 
spin (iTchekhovskov et akl 1201 ll) . instabilities at the jet-disk 
interface limit how high rjet can be in practice. Thus, even 
for black hole spin-powered jets, Ljet may be used to estimate 
M, and thus address the existence of an event horizon in M87*. 

In Section |2] we describe how the jet power and M are 
related, and obtain an estimate for the latter in M87*. The size 
constraints placed by mm-VLBI are summarized in Section 
[3 and the relevant observational flux limits are presented in 
Section H) The associated constraints upon the existence of 
a photosphere are discussed in Section |5] Our conclusions are 
summarized in Section|6l 


2. ESTIMATES OF M 


2.1. Relating M to Jet Power 

All current models for launching relativistic jets require an 
accretion disk. In the simplest models, the jet is merely 
the innermost relativistic part of a magnetocentrifugal wind 
flowing out from the surf ace of the disk. In such models (e.g., 
IBlandford & Pavn'^ll982^ ■ the jet and the lower-velocity wind 
are ultimately powered by the gravitational potential energy 
released by the accreting gas. Thus the jet luminosity is limited 
by the overall energy efficiency of the disk. 


Tjet Tjet Tfadiation 


( 1 ) 


where the disk efficiency lydisk is determined by the radius 
of the inner edge of the disk, i.e., the radius osco of 
the innermost stable circ ular orbit (ISCO) of the space-time 
dNovikov & Thorn3ll973l) . 


?7disk = 1 - t/1 - —• (2) 

V 3risco 

The efficiency pdisk varies from 0.057 for a non-spinning black 
hole up to 0.42 for a maximally spinning black hole. A typical 
value is probably ~ 0.1-0.2. This means that, for disk- 
powered jets, the mass accretion rate implied by a given jet 
luminosity is roughly 

M«10%. (3) 

Alternatively, the jet may be powered by black hole rotation 
(IBlandford & Znai^ll977h . In this case, the jet is launched by 
large-scale ordered magnetic fields that penetrate the horizon, 
and the total jet power is 


Liet = — 

J 47rc " 


(4) 


where VIh = a^:cl2r+ is the angular velocity of the horizon, 
located at r+ = {GM/c^)(\ + yjl-al) where a* = a/M is 
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the dimensionless spin of the black hole, $ is the magnetic 
flux threading the horizon, and k « 0.045 is a dimensionless 
coefficient that varies mo destly with black hole spin (see 
iTchekhovskov et ^120101 for a numerical calculation of this 
dependence). As Eq. dUl shows, the jet power increases 
quadratically with both black hole spin and magnetic flux. The 
former is limited by the condition a, < 1, and the latter is 
limited by the requirement that the magnetic held through the 
horizon has to be confined by the accretion flow. As a rough 
estimate, one could say that the accreting gas is virialized and 
exerts a ram pressure of order 

T’acc ~ pvl , (5) 

where is the local Keplerian velocity. At the disk-funnel 
interface, assumed to occur at a cylindrical radius rgq, the ram 
pressure must balance the magnetic pressure within the jet, and 
thus 

2 \2 2kcTr^ 4 2 

^jet ~ V^^^eq^eqj ~ ^eqPeq^fc,eq i (^) 

where Bgq, Peq, and v^^gq are the magnetic field strength, gas 
density and Keplerian velocity at rgq and we have assumed 
a* ~ 1. 

At the disk-funnel interface, the heavier gas is supported 
on top of an otherwise buoyant magnetic held, and thus a 
balanced configuration is generally unstable. The interchange 
instability, a close relative to the more commonly discussed 
Raleigh-Taylor instability, results in the growth of gas Angers 
that interpenetrate the magnetosphere, allowing accretion to 
occur on time scales comparable to the Keplerian period at 
the interface (IS pruit et aki 119951 : llgumenshchev et 120031 : 

iLi & NaravanI l2004li . Therefore, even in the presence of a 
strong funnel field, the typical accretion rate is 

M « 47rr“qPgqVi:,gq , (7) 

differing from the previous estimate by only a factor of order 
unity. As a result, 

M«10%, (8) 

i.e., this approximate calculation gives an estimate for M not 

v ery different from the case of a dis k-powered jet (Eq. [3]l. 

iGhosh & Abramowiczl(ll997l) and lLivio et al.l(ll999h studied 
in greater detail the balance between an accretion disk and 
a magnetic flux bundle confined around a central black hole 
and concluded that, for the case of a standard radiatively 
efficient thin accretion disk, the likely jet efficiency is 
significantly smaller than 10%, which means that M is likely 
to be substantially larger than lOLjgt/c^. The situation is, 
however, different in the case of systems like M87 which have 
hot advection-dominat ed accretion flows (for a review, see 
lYuan & Naravanll2014^ . 

Rec ent GRMHD simu la tions (e. g.. ITchekhovskov et akl 
201 ll iMc Kinn ey et al.l 120121 iNaravan et al.l 120121 
Sadowski et akl 120131) have shown that substantial magnetic 
flux can be confined around a black hole by a hot accretion 
flow. The field strength is larger than simple estimates suggest 
partly because of geometrical factors related to the shape 
and dynamics of the accretion flow and partly because of 
relativistic correction s in the vicinity of the black hole horizon 
(It chekhovskovll20T5h . The net result is that the jet efficiency 
can be up to a factor of 10 larger than naive estimates suggest. 
Indeed efficiencies exceeding 100% are possible if the black 


hole spins at close to the maximal rate and the magnetic 
flux achieves its maximum value via a magnetically arrested 
disk (MAD, INaravan et akl 120031: ITchekhovskov et d11201 ll) 
configuration. Allowing for these effects, we thus have 


where the value we give for the coefficient is highly conserva¬ 
tive. In the following, we make the conservative assumption 
thatM = Ljet/2c^. 

2.2. Estimates of Jet Power in M87 

A variety of estimates of M87’s jet power may be found in 
the literature, covering a wide variety of distances from the 
black hole and thus timescales probed. All the estimates are 
consistent with a jet power of roughly 10"*^ ergs”'. 

Surrounding M87 is an ext e nded radio-bright structur e 
(iBolton et akl [1941 ll^^ [1951 iBaade & Minkowskil [T954I) . 
reaching radial distances of ne arly 30 kpc, and be lieved to be 
powered by the central AGN (lOwen et akl 120001) . Estimates 
of the jet power may be obtained by estimating the energy 
required to grow the radio halo and dividing it by the buoyancy 
times cale; this gives Ljet « few x 10"'"' ergs”' (lOwen et alJ 
l2000h . Recent efforts to estimate the synchrotron age of 
the halo using LOFAR observations give roughly 40 Myr, 
resulting again in Ljet ~ 6-1 0 x 10"'"' ergs”', depen ding on 
the assumed particle content (Ide Gasnerin et akl 120121) . These 
necessarily represent estimates averaged over the buoyancy 
time at 10 kpc, roughly 10 Myr. As a result, inferring 
the current instantaneous jet power requires some assumption 
regarding the recent history of activity in M87. 

Like many radio-loud AGN, on kpc scales M87 exhibits X- 
ray cavities and regions of enhanced emission associated with 
shocks. These are presumed to be driven by the AGN jet, and 
thus provide an additional measure of jet power. Estimates 
of the shock energetics alone require a power source of 2.4 x 
10"'^ ergs”' (iForman et akll200-'^ 120071) . The power inferred 
from the nearby cavity inflation is sensitively dependent on the 
gas pressure profile, and estimates range fr o m 10"'^-10"'" ' ergs” ' 
(I Young et akl 120021 : iRaffertv et alJ l2006t lAllen et ^ 120061) . 
Again, these are time-averaged over the buoyancy time, in this 
case roughly ~ 1 Myr because of the smaller scale (1 kpc), with 
the attendant caveats regarding variability. 

Knot A is a bright optical feature within the jet at 
approximately 0.9 kpc. While comparable to the distance 
of the closest cavities, knot A is one of a number of 
superluminal features within the jet, with apparent velocities of 
up to 1.6c, implying that it moves relativistic ally (iMever et akl 
120131) . Interpreting knot A as an oblique shock within the 
jet results in a jet power es timate of (1-3) x 10"'"' ergs”' 
(iBicknell & BegelmanI [l99^; values much in excess of 
this are expected to over-produce the emission at knot A 
([Reynolds et aljn~996h . Equally importantly, the time delays 
associated with estimates from knot A are roughly 2x10’' yr, 
three orders of magnitude shorter than those associated with 
the large-scale radio and X-ray morphology. 

Finally, located 60 pc from M87* is the HST-1 complex, 
comprised of stationary and s uperluminal components with ap¬ 
paren t velocities of up to 6c (iBiretta et aklll99^ iGiroletti et akl 
l2012b . At these distances, HST-1 provides the most contempo¬ 
raneous estimates of the jet power, with a time delay of roughly 
30 yr. If identified with a recollimation shock, the station¬ 
ary component of HST-1 implies a jet power of w 10"'"' ergs”' 
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(IStawarz et al.ll200^ . This is grossly consistent with estimates 
obtained when the shock stru cture is modeled in more detail 
(iBromberg & Levinsonll200^ . 

Overall, the typical estimate for the bolometric power 
of M87’s jet is ^ 10"^'* ergs“', with perhaps a factor of 
few uncertainty. The implied mass accretion rate is M ^ 
IQ-^Moyr"' 

3. 1.3MMVLBISIZEOFM87 

The detection of emission at the core of M87 on the 
scale of several times the Schwarzschild radius is robust. 
Observations on an array including telescopes in California 
(CARMA), Hawaii (JCMT) and Arizona (SMT) provided 
projected baseline lengths ranging from 6 x 10^Ato3.4x lO'^A. 
Eirm detection of M87 on the longest baselines alone requires 
the existence of compact structure with an upper size limit 
of ~ 60 microarcseconds, corresponding to a diameter of ^ 
16GM/c^ at the presumed mass and distance of M87. When 
modeled as a single c ircular Gaussian brightness distribution 
(iDoeleman et al.l[2012h . the Eull Width Half Maximum is 11 ± 
XGMjc^ with calibration errors included in the S-cr parameter 
estimate. 

Other instrumental factors can affect the estimated 1.3 mm 
VLBI size but were determined to be negligible. These 
include the stability of Hydrogen Maser frequency standards 
used at each VLBI station and the polarization purity of the 
telescope receivers which were conhgured for Left Circular 
Polarization. On short (si 1 s) timescales, the Masers 
were compared to ultra-stable quartz oscillators, and longer- 
term stability was determined through timing comparisons 
with GPS. Polarization characteristics at each antenna were 
determined by injecting known polarization signals into the 
receiver feeds at the single dish sites, and through separate 
polarization calibration observations at CARMA. 

The main source of uncertainty in the true size and shape of 
the compact 1.3mm emission is the limited VLBI baseline cov¬ 
erage, which is insufficient to exactly determine the brightness 
morphology. The VLBI data can also be fitted, for example, 
with a uniform disk of diameter 16.4GM/c^, comparable to 
that inferred from the the long-baseline correlated flux mea¬ 
surements alone. In this work we adopt the circular Gaussian 
diameter of WGMjc^ with the proviso that future EHT obser¬ 
vations with improved Global coverage will be able to image 
and model M87 in detail. The disk and more complex mod¬ 
els typically predict “nulls” in the VLBI signal as a function of 
baseline length that are not yet seen, but which may emerge as 
E HT sites enabling longe r baselines are added to the array. 

IDoeleman et al.l (1201 2b associate the 1.3mm compact emis¬ 
sion with the ISCO of the M87 black hole, enlarged by the 
strong gravitational lensing that occurs at small radii. Loca¬ 
tion of this emission close to, or at, the black hole is sup¬ 
ported by phase referenced multi-frequency VLBI observations 
at longer wavelengths, which show the core of M87 to shift to- 
wards a converge nt point with shorter observing wavelength 
jHada et al.ll201 ll) . Extrapolation of this shift to the observed 
size at 1.3mm places the emission within a few Schwarzschild 
radii of the black hole, making the derived 1.3mm size a good 
estimate for the photosphere of a putative surface. 

While it is possible that the 1.3mm emission arises some 
distance from the black hole, the EHT-derived size can still 
serve as an upper limit on the spatial extent of an intrin¬ 
sic photosphere. Theoretical models of electromagnetic jets. 


both semi-analytical estimates and numerical simulations, uni¬ 
formly show an expanding jet width prohle with increased dis¬ 
tance from the central engine. These expectations are quan¬ 
titatively consistent wit h jet observations on scales covering 
six orders of magnitude (iHada et alJl2013tlNakamura & Asadal 
120131) confirm a parabolic jet shape to within ^ IQGMjc^ of 
the black hole. There is, in short, no compelling mechanism 
to suggest that an electromagnetic jet launched from a photo¬ 
sphere later converges to a markedly smaller size. 

4. THE SED OE THE NUCLEAR AGN SOURCE IN M87 

Tables [T] and |2] list the optical and near-infrared SED 
data that are used in Ligur es [T] and Table p| lists the 
measurements made in the iPerlman et al.l (1201 ih five-year 
campaign to measure nuclear variability. Table [T] lists the rest 
of the data. 

Most flux measurements were made with the Hubble Space 
Telescope (HST) and various imaging cameras. Also, unless 
otherwise noted, nuclear flux measurements were made using 
a correction for the galaxy light under the central PSL that is 
base d on measuring and fitting the profile at larger radii (see, 
e.g.. iLauer et alJ[T992t iSnarks et al.lll996^. V arious prohle ht 
models s uch as Sersic funct ions (ISersiclll96§) and ones with a 
core (see ILauer et £0119951) are used. Also unless otherwise 
noted, the innermost jet knot HST-1 is not included in the 
measurements, although fai nt extensions of the jet at smaller 
radii are included (see, e.g.. iSnarks et al.lFl996l: l^rlman et al.l 
120011) . The measurements are discussed in their source papers. 
A few details are noteworthy here: 

An estim ate byl^ung et alJ(ll978h of V ~ 16.69 ±0.05 (see 
ILauer et ^1199^ corresponding to Tv = 0.80 ±0.04 mJy is 
consistent with the results in Table [T] but is not used in this 
paper. The epoch was 1975 Mar 14 - 1977 May 22 and the 
Palomar Observatory 200-inch Hale telescope was used. We 
omit this measurement because the seeing LWHM was 1" to 
2" - hence the exclusion of jet knots is uncertain - and because 
the galaxy subtraction was not consistent with our HST-era 
understanding of the c uspy cores of elliptical galaxies (e.g., 
ILauer et aklfl^flO^. 

The ILauer etal.1 (119921) HST PC P785LP images (Aeff = 
8900 A) were taken with the unrefurbushed HST but allow a 
clear separation of the nuclear AGN from the innermost jet 
knots N1 and M at ~07l and ^ 0"18 from the center (see their 
Ligure 5). These knots are includ ed in most HST ph otometry 
and all ground-based photometry. ILauer et aP (119921) measure 
them separately and show that, at 8900 A and on 1999 Jun 1, 
the nucleus contributed 89% of the combined flux including N1 
and M. We have used their measurement of the nuclear flux at 
8900 A. They also measure the total flux of the nucleus plus 
N1 plus M in the L555W band at almost the same time as the 
L785LP measurements. We assume that the fraction of this 
total that comes from the nucleus is also 89% and correct their 
total measurement to a nuclear estimate of 0.637 mJy. This 
correction has not been applied to any other measurements, 
because we cannot know the variability of the nucleus, ofNl, 
and ofM separately. All these other measurements include the 
innermost part of the jet interior to knot HST-1. 

We repeated the L55W galaxy subtraction to get the 
estima ted error ±0.059 on t he L555W flux. 

The ISnarks et al.l (119961) measurements were made with 
pre-COSTAR HST using the ground-based photometry by 
lYoung et al.l (119781) as the galaxy model. The latter is probably 
not a problem; they note that corrections for galaxy light are 
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Table 1 

Ultravoilet-Infrai'ed SED of the M87 Nuclear Source 


A 

(fim) 

(1) 

logi/ 

(Hz) 

(2) 

F^ 

(mly) 

(3) 

CTF 

(mly) 

(4) 

Epoch 

(5) 

Program 
HST ID 
(6) 

Camera 

(7) 

Source 

(8) 

0.1255 

15.38 

0.162 

0.032 

1991 Apr 6 

1228 

FOC 

et al. (1996) 

0.1507 

15.30 

0.081 

0.016 

1991 Apr 5 

1228 

FOC 

et al. (1996) 

0.1585 

15.28 

0.068 

0.014 

1991 lun 23 

1517 

FOC 

^gaifcetal. (19J)()) 

0.23 

15.11 

0.158 

0.032 

1991 Apr 5 

1228 

FOC 

^g^x et al. (1996) 

0.25 


0.1153 

0.0014 

2003 Mar-31 

9454 

ACS HRC 

Maoz et al. G005) 

0.25 


0.1119 

0.0014 

2003 Dec 10 

9454 

ACS HRC 


0.2475 


0.089 

0.009 

1999 May 17 

8140 

STIS NUV-MAMA 

Chiaberue et al. (2002) 

0.299 


0.198 

0.005 

1998 Feb 25 

6775 

WFPC2 PC 

This Daoer: Perlman et al. (2001) 

0.33 


0.1759 

0.0019 

2003 Mar 31 

9454 

ACS HRC 

Maoz et al. (2005) 

0.33 


0.1743 

0.0019 

2003 Dec 10 

9454 

ACS HRC 


0.3708 

14.91 

0.457 

0.091 

1991 lun 23 

1517 

FOC 

^g^x et al. (1996) 

0.456 


0.493 

0.012 

1998 Feb 25 

6775 

WFPC2 PC 

This paper’ Perlman et al. (2001) 

0.5017 

14.78 

1.02 

0.20 

1991 Apr 5 

1228 

FOC 

Jlg^x et al. (199()) 

0.541 


0.637 

0.059 

1991 Feb 24 

1105 

WF/PC PC 

Lauer et al. (1992): this paper 

0.600 


0.787 

0.020 

1998 Feb 25 

6775 

WFPC2 PC 

This paper’ Perlman et al. (2001) 

0.801 


1.063 

0.027 

1998 Feb 25 

6775 

WFPC2 PC 

This nane^^erlman et al. (5651) 

0.801 


0.759 

0.076 

1999 May 11 

8140 

WFPC2 

Chiaberse et al. (1999^002) 

0.8900 


0.81 


1991 lun 1 

3242 

WF/PC PC 

Lauer et al. ^92) 

1.2 


2.45 

0.69 

1993 lun 1 


ESO/MPI 2 m 


1.60 


3.02 

0.15 

1997 Nov 20 

7171 

NICMOS NIC2 

Baldi et al. (2010) 

2.201 


2.28 

0.60 

1994 Apr 3 


UKIRT 3.8 m 

Snarks et al. (1996) 


Note. — NOTES - Columns (1) and (2) are the effective, pivot, or median wavelength and frequency of the photometric bandpass. 
Columns (3) and (4) are the measured flux and its l-cr error. Because the nuclear source flux is time-variable. Column (5) gives the 
epoch of the observation. This is not used in the present paper: all SED points are plotted in Eigure 1 to provide an illustration of the 
amplitude of variation. This amplitude is in general larger than the errors of the SED measurements. Column (6) provides the HST 
proposal ID, and Column (7) lists the HST camera or the telescope used. Column (8) gives the source of the measurement. 


Table 2 _ _ 

Ultravoilet-Optical SED of the M87 Nuclear Source from IPerlman et al.l 1201 Ih 


Fv (mJy) 
F606W 
0.589 /rm 
(1) 

Fv (mJy) 
F330W 
0.336 /rm 
(2) 

Fv (mJy) 
F250W 
0.272 /rm 
(3) 

Fv (mly) 
F220W 
0.226 ^m 

(4) 

Epoch 

(5) 

Program 
HST ID 

(6) 

Camera 

(V) 

0.671 ±0.007 

0.305 ±0.002 


0.146 ±0.011 

2002 Dec 7 

9705 

ACS HRC 

0.630 ±0.006 

0.305 ±0.002 


0.146 ±0.011 

2002 Dec 10 

9493 

ACS HRC 

0.478 ±0.005 



0.137±0.011 

2003 Nov 29 

9829 

ACS HRC 

1.066 ±0.011 

0.475 ±0.002 

0.306±0.014 

0.226 ±0.014 

2004 Nov 28 

10133 

ACS HRC 

1.306±0.013 


0.363 ±0.010 


2004 Dec 26 

10133 

ACS HRC 

0.891 ±0.009 


0.280 ±0.009 


2005 Feb 9 

10133 

ACS HRC 

1.037 ±0.010 


0.328 ±0.010 


2005 Mar 27 

10133 

ACS HRC 

0.932 ±0.009 

0.446 ±0.003 

0.274 ±0.009 

0.217 ±0.009 

2005 May 9 

10133 

ACS HRC 

0.839 ±0.008 


0.273 ±0.009 


2005lun 22 

10133 

ACS HRC 

0.639 ±0.006 


0.192 ±0.007 


2005 Aug 1 

10133 

ACS HRC 

0.735 ±0.007 

0.349 ±0.002 

0.217 ±0.008 

0.170 ±0.008 

2005 Nov 29 

10133 

ACS HRC 

0.756 ±0.008 


0.234 ±0.008 


2005 Dec 26 

10617 

ACS HRC 

0.631 ±0.006 


0.201 ±0.008 

0.160 ±0.008 

2006 Feb 8 

10617 

ACS HRC 

0.780 ±0.008 


0.232 ±0.008 


2006 Mar- 30 

10617 

ACS HRC 

0.862 ±0.009 

0.372 ±0.002 

0.193 ±0.007 

0.1564 ±0.007 

2006 May 23 

10617 

ACS HRC 

1.370±0.014 

0.636 ±0.003 

0.323 ±0.009 


2006 Nov 28 

10910 

ACS HRC 

1.006 ±0.010 


0.276 ±0.009 


2006 Dec 30 

10910 

ACS HRC 

1.292 ±0.017 




2007 Nov 25 

11216 

WFPC2 WFC 


Note. — Notes - The first four columns list fluxes (mJy) for the four filters and their pivot wavelengths as 
listed in the table headers. Columns (5) - (7) are explained in their headers. 


small. If the former is a problem, then light at large radii in 
the PSF could be missed and the measurements could be lower 
limits on fluxes. However, we emphasize that this problem 
is almost certainly much smaller than the intrinsic variability 
a mplitude of the nucle ar source. 

ISparks et al.l (119961) measured the nuclear AGN in J and 
K bands using the ESO/MPI 2 m telescope and the United 
Kingdom Infrared Telescope (UKIRT: 3.8 m), respectively. 


Galaxy light underlying the AGN was removed using an 
iterative procedure to separate the stellar brightness profile 
from the central source. That source was measured inside a 
2"8-diameter aperture, i.e., one that is large enough to include 
jet knot HST-1. The knot is not visi ble in their image , 
and the higher-resolution HST images in ISparks et all (Il996h 
suggest that it contributed negligibly at that time. There are jet 
components closer to the center ('e.g.. lLauer et al.l^l992^ : these 











































6 


BRODERICK ET AL. 


are inc luded even in most HST measurements. We therefore 
use th elSnarks et al.l (^1996^ measure ments as published. 

The lChiaberge et al.l (Il999ll2002h measurements were made 
by subtracting from the nuclear flux the galaxy flux at radius 
0"23 in the R band and at radius 0"17 in the ultraviolet band. 
Since the g alaxy’s profile is a shallow power law inside the 
cuspy core (iLauer et al.lll992h . the resulting fluxes are strictly 
speaking upper limits. However, the galaxy contribution to the 
measurements is smaller than the amplitude of the variability 
of the M 87 nuclear source. 

The four HST WEPC2 points from t his paper were measu red 
by E.S.R following the procedures of iPerlman et al.l (1200 Ih in 
Eebruary of 1998. Knot HST-1 was about 1% of the central 
source at this time and is not included in the measurement. 
Prior to measuring the nuclear flux, a galaxy model (made with 
the IRAE tasks ellips e and bmodel) was s ubtracted from 
each of these images. The lPerlman et al.f(l201 Ih measurements 
listed in Table|2]were made in the same way. 

We have n ot included the UV photometry of the nucleus 
published by iMadridI (l2009h because all of those fluxes have 
been referenced to a single wavelength based on an unpub¬ 
lished sp ectral model. Instead , we use the measurements pub¬ 
lished bv iPerlman et alJ (1201 Ih . which incl ude all of the ep ochs 
between 2004-2006 that were included in IMadridI (1200^ but 
does not apply any spectral model in their calibration. I t should 
be not ed that when the spectral referenci ng applied by MadridI 
(l2009f) is removed , the measurements of lPerlrna n et al. T201 ll) 
and lMad^ (l2009h agree to within the errors (lMadridll201 Ih . 

5. PHOTOSPHERE EMISSION 

Were a photosphere present in M87* instead of a horizon, it 
would be heated by the constant deposition of kinetic energy 
from the accretion, resulting in an additional component in 
the spectrum. Since the infall time is only logarithmically 
dependent upon the photosphere redshift, it is natural to assume 
that the surface will reach a steady state in which the emitted 
luminosity and the impinging kinetic power from the accretion 
flow are balanced. 

The gravitational binding energy liberated per unit time by 
the infalling material necessarily depends on the size of the 
putative photosphere, approaching Mc^ for radii comparable to 
that of the horizon. That is, typically, 

GM . , 

( 10 ) 

where R is the radius of the surfaceQ 

The radius of any such surface is strongly limited by the 
1.3mm VLBI observations discussed in Section [3] These 
already constrain the 1.3mm emission region to within the 
vicinity of the photon orbit, expected to have a diameter 
of roughly 10.4GM/c^ in projection and depending very 
weakly on the assumed black hole spin0 Moreover, at 
millimeter wavelengths and longer M87 still exhibits the flat 
radio spectrum indicative of a self-absorbed synchrotron jet 

^ While quantitative variations from the Newtonian expressions do occur 
for high-redshift objects, these are qualitatively similar, differing by a factor of 
order uni ty that depends primar ily on the dynamical state of the infalling gas 
(see, e.g.. (Broderick et al.l2009l and Equation[3. 

^ While it is true that the projected photon orbit viewed from the equatorial 
plane can shrink to 9GM/c^ for maximally spinning black holes, this is not the 
cas e at small inclinations; at t he estimated inclination of < 25° of the M87*’s 
jet 4Heinz & Begelm^ll997f) . and thus presumably spin inclination, spin has 
negligible impact on the projected photon orbit size. 


(see, e.g., iBlandford & Koni^ Il979ll . This interpretation 
receives strong support from the evolution of the radio core 
location with wavelength, with the offset evolving o c ^0-94ibo.09^ 
presu mably due to the shrinking of the photosphere (iHada et al.l 
1201 ll) . Thus, at infrared and optical wavelengths the size 
of the photosphere is expected to be smaller than the ~ 
llGM/c^ limit from the mm-VLBI observations described in 
Section [3 - if the optically thick to thin transition completes 
at a wavelength of 1.1mm or shorter, this will lie within the 
projected photon orbit. This conclusion is only strengthened if 
the emission arises from a significant distance from the central 
engine, as described in Section[3] 

However, the efficiency with which the jet is driven depends 
on the surface in a similar fashion, becoming larger for more 
compact objects. Eor disk driven jets, the surface and jet 
efficiences are approximately equal (see Equation[Tli, implying 
that Lsurf ~ Tjet. Eor black hole rotation driven jets the 
efficiency scales roughly as due to the competition between 
the lower H// and higher 4> as the surface increases in size (see 
the Appendix). As a result, once again Lsurf ~ Tjet. with an 
additional factor of order unity that depends only very weakly 
on the photosphere size. Thus, here we make the conservative 
assumption that Lsurf ~ Tjet/2 as described in Equation|9] 

Since the size constraints placed by EHT observations at 
1.3mm already limit the size of any putative photospheric 
emission in M87* to scales comparable to, or smaller than, the 
photon orbit, the photosphere will necessarily be thermalized 
by the strong lensing, resulting in a gravitational analog 
of the standard blackbody cavi ty (see the discussion in 
iBroderick & NaravanI 120061120071) . Thus, surface should be 
temperature as seen by a distant observer is related to the 
accretion rate via 


4nRlaT^ ^ ^ 


\%TiRla) 


1/4 


( 11 ) 


where a is the Stefan-Boltzmann constant and Ra « 5GM/c^ is 
the apparent radius of any object lying within the photon orbit. 
Note that the temperature estimate depends only weakly on 
M, and thus on the uncertainties associated with the jet power 
and jet efficiency. Eor the accretion rate inferred in Section |2] 
for M87*, Lx) « 8.5 x 10^ K, which means that the surface 
emission will peak in the optical band. The observed spectral 
flux from the surface is given by 


7v,ph 


R^ 2hv^/c^ 
Jf- ehvIkT^ _ 1 ’ 


( 12 ) 


where D is the distance to M87. This predicted flux may be 
directly compared to observations and thereby we can hope to 
constrain whether or not M87* has a surface. 


6. DISCUSSION 

A direct comparison of the putative photosphere spectrum 
to the measured fluxes listed in Table [T] is shown in Eigure 
[U Because of the intrinsic uncertainty in M, we show 
predicted spectra spanning a wide range of M, corresponding 
to variations in the jet efficiency or measured jet power. In all 
cases, photospheric emission is excluded at high confidence. 
This conclusion is further supported by Eigure |2] which 
shows the ratio of the measured fluxes to that implied by a 
photosphere as a function of M. At all physically reasonable 
accretion rates the measured fluxes lie at least an order of 
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1 0.1 



Figure 1. Inferred infrared-optical spectrum from a putative photosphere for 
f-photo = 10^^ ergs^*, 10^** ergs^*, and 10^^ ergs“* (solid lines, from bottom to 
top) compared to the empirical limits on the SED of M87* listed in Tables [T] 
and |2] (colored dots). The dark-red dotted line indicates the estimated limit 
arising from the intrinsic jet spectrum, and thus presumably the confusion 
limit. All points in Table |2| are plotted in order to provide some indication of 
the variability of the source. Those data that provide the constraints presented 
in Figure |2] are shown in red. Note that photospheric emission at a level of 
10'*“* ergs^*, as expected based on current estimates of the jet power (§2), is 
ruled out with high confidence. 

magnitude belo'w that required in the absence of an event 
horizon. 

Improvements in the measured flux limits 'will necessarily 
result in increased confidence with which any photospheric 
emission may be ruled out. However, the intrinsic emission 
from the observed jet component provides a natural limit. 
Assuming a spectral index of 1.2, consistent 'with a transition 
from optically thick to optically thin near 1mm, and the 
optical flux limits, we estimate the corresponding floors on the 
infrared-optical spectrum and flux ratio, sho'wn by the dotted, 
dark-red lines in figures [T] and |2l belo'w which the emission 
from the jet-launching region must be modeled in detail. Note 
that this is rather uncertain; variations in the spectral index 
bet'ween the core and larger scale jet emission arising, e.g., 
from an evolution in the relativistic particle distributions, could 
easily harden the core spectrum sufficiently to fully account for 
the vast majority of the fluxes shown in Figure[T] Nevertheless, 
even if this is not the case, it appears that the current flux limits 
are 'within an order of magnitude of the confusion limit, below 
which careful modeling of the jet emission 'will be required to 
obtain stronger constraints. 

Despite the astrophysical uncertainty in the relationship 
bet'ween the jet po'wer and mass accretion rate, we have sho'wn 
in this paper that, within the context of current jet launching 
paradigms, 'we can rule out the existence of an observable 
photosphere in M87* 'within which the kinetic energy of the 
accreting gas is deposited. The implication is that the kinetic 
energy of the gas is advected past an event horizon, beyond 
which it is no longer visible to distant observers. In other 
words, M87* must have an event horizon. 



Figure 2. Ratio of the measured fluxes to that inferred from a putative 
photosphere. The thick black lines indicate the constraint implied by the 
individual flux measuremerrts shown in red in Fig[T] corresponding to those 
that dominate the constraint at some M in the range presented. The line widths 
indicate the la upper limit and the grayed region is the area excluded by the 
aggregate. The dark-red dotted line indicates the estimated limit arising from 
the intrinsic jet spectrum, and thus presumably the confusion limit. 

This argument may be evaded in a variety of 'ways, though 
all require speculative new physics and most invoke black hole 
alternatives. First, M87* could fail to reach steady state, pre¬ 
venting the use of the accretion rate as a proxy for the sur¬ 
face luminosity. This requires un physically high heat cap aci- 
ties (which the gravastar model of iMazur & Mottolall2004l ap¬ 
parently doe s possess, seelC haplinel 2005[) . already excluded in 
many cases (iBroderick & Naravanl2007h. or exotic alternatives 
like s uitably tuned wormholes (e.g.. iDamour & SolodukhinI 
I2007h . Second, the surface could fail to thermalize, though 
this is strongly argued against by the compact emission ob¬ 
served at millimeter wavelengths, and may be excluded alto¬ 
gether by the development of a baryonic atmosphere. Third, 
the efficiency factor could be much larger than implied by cur¬ 
rent jet launching models, requiring mechanisms qualitatively 
different from those currently favored. For example, if M87* 
were a rapidly spinning compact object 'with a strong magnetic 
field, the accreting gas could be stopped at a magnetospheric 
radius and flung out in a jet 'without any gas reaching the un- 
derlying photosphere (e.g., the neutron star propeller model of 
llllarionov & Sunvaevlll975]) . In any such model, the magne¬ 
tospheric radius can be only slightly larger than the surface 
of the central object (because of the VLBI constraints on the 
angular size of the source), and it seems likely that a lot of ac¬ 
creting gas 'will penetrate through the magneto sphere and reach 
the surface (e.g.. lKulkarni & Romanoval2008h . In addition, the 
material that is stopped and flung out by the magnetosphere is 
likely to dissipate a good fraction of its kinetic energy 'where 
it meets the rotating magnetosphere, resulting in thermalized 
radiation not very different from that expected from a photo¬ 
sphere. We regard this propeller model as strictly physically 
possible but astrophysically implausible. Thus, our argument 
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will be strengthened in the near future as key components to 
the current jet paradigm are critically tested by mm-VLBI ob¬ 
servations of M87. 


Based on observations made with the Event Horizon 
Telescope and the NASA/ESA Hubble Space Telescope. 
The Event Horizon Telescope is supported through grants 
from the US National Science Eoundation, by the Gordon 
and Betty Moore Eoundation (GBME3561), and through 
generous equipment donations from the Xilinx Corporation. 
The NASA/ESA Hubble Space Telescope observations were 
obtained at the Space Telescope Science Institute, which is 
operated by AURA, Inc., under NASA contract NAS 5-26555. 
These observations are associated with program numbers 1105, 
1228, 1517, 3242, 6775, 7171, 8140, 9454, 9493, 9705, 
9829, 10133, 10617, 10910, and 11216. A.E.B. receives 
financial support from Perimeter Institute for Theoretical 
Physics and the Natural Sciences and Engineering Research 
Council of Canada through a Discovery Grant. RN received 
partial supported from the NSE via grant AST1312651 
and NASA via grant NNX14AB47G. RN also thanks the 
Perimeter Institute for hospitality while some of this work 
was carried out. The authors are grateful to Sera Markoff 
and Brian McNamara for helpful conversations. This work 
would not have been practical without extensive use of 
the NASA/IPAC Extragalactic Database (NED), which is 
operated by the Jet Propulsion Laboratory and the California 
Institute of Technology under contract with NASA. We also 
made extensive use of NASA’s Astrophysics Data System 
bibliographic services. JK’s work was supported in Texas 
by the Curtis T. Vaughan, Jr. Centennial Chair in Astronomy 
and in Germany by a Eaculty Research Assignment from 
the University of Texas and by the Max-Planck-Institut 
fiir Extraterrestrsiche Physik (MPE), Garching-by-Munich, 
Germany. JK warmly thanks Director Ralf Bender and the 
staffs of the MPE and the Universitats-Sternwarte, Ludwig- 
Maximilians-Universitat, Munich for their hospitality and 
support during the 2014 visit when most of his work on this 
paper was done. 

Eacilities: EHT; James Clerk Maxwell Telescope, Arizona 
Radio Observatory Submillimeter Telescope, Combined Array 
for Research in mm-wave Astronomy, the Submillimeter Array 
HST: EOC, HST: WEPCl, HST: WEPC2, HST: NICMOS, 
HST: ACS, HST: STIS, UKIRT 3.8 m telescope, ESO/MPI 2 
m telescope. 


APPENDIX 

A. BLANDFORD-ZNAJEK JET EFFICIENCY FOR SURFACES 

The relationships between accretion rate and both surface 
luminosity and jet power are dependent upon the assumed 
surface size. In detail this depends upon the specific jet 
launching mechanism assumed. Here we consider the question 
in the context of Blandford-Znajek jets, for which the jet 
luminosity is related to the properties of the central object 
via Equation |4] Eor an object with a large surface there are 
two important modifications: the angular velocity is limited to 
Hh < c/R, and <1> < Thus, 

Ljet < « 27rk^Mc^ (Al) 


where we assumed M « Pea}’k,eq and equipartion near the 
surface (i.e., req « R). This must be compared to 

GM ■ . ■ 0 

Tsuif ~ ^ ^^k,eq > 

which gives 

Tsurf^, 1 Vfc.eq ^ 

Ljet 27rk c \ c^R J 

This is a weak function of R. Eor the range of surface radii 
consistent with the recent mm-VLBI limits, it introduces at 
most a factor of two. 


(A2) 


(A3) 
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